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IDENTIFICATION OF LUMINOUS INFRARED GALAXIES AT 1 < Z < 2.5 
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ABSTRACT 

We present preliminary results on 24 /im detections of luminous infrared galaxies at z > 1 with the 
Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole 
and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array 
Camera (IRAC) between 3 and 9 /im. The positional accuracy of < 2" for most MIPS/IRAC detections 
provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from 
the Deep Extragalactic Evolutionary Probe survey, we identify 24 /im sources at z > 1 in the EGS, while 
the combination of the MIPS/IRAC observations with BVRIJHK ancillary data in the Lockman Hole 
also shows very clear cases of galaxies with photometric redshifts at 1 < z < 2.5. 

The observed 24 /im fluxes indicate infrared luminosities greater than 10 11 L Q , while the data at shorter 
wavelengths reveal rather red and probably massive (A4 > A4*) galaxy counterparts. It is the first time 
that this population of luminous objects is detected up to z ~ 2.5 in the infrared. Our work demonstrates 
the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how 
the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic 
history not explored by previous infrared experiments. 

Subject headings: galaxies: high-redshift — infrared: galaxies — cosmology: observations 



1. INTRODUCTION 

In the past few years, deep observations in the infrared 
and submillimctcr revealed a population of high redshift 
galaxies emitting the bulk of their luminosity between 8 
and 1000 /im, whose cosmological significance had been 
previously missed or sevcrly underestimated by optical 
surveys because of extinction effects due to dust. These 
sources, likely the analogs of the local Luminous and 
Ultra-Luminous InfraRed Galaxies (respectively LIRGs: 
10 n L o < L IR =L[8 - 1000/im] < 1O 12 L , and ULIRGs: 
Lir > 1O 12 L ) played a crucial role through the cosmic 
ages. At mid-infrared wavelengths, results from the In- 
frared Space Observatory (ISO) mission show that the 
evolution of these powerful objects noticeably contributed 
to the global star formation history at < z < 1 (Elbaz 
et al. 1999; Flores et al. 1999). Observations with SCUBA 
at 850 /im also led to the discovery of extremely bright in- 



frared sources mainly located at z > 2 and characterized 
by a space density 1000 times larger than observed in the 
local Universe (e.g., Chapman et al. 2003). 

From the purely observational point of view though, 
very little is currently known about the potential impor- 
tance of dusty sources at 1 < z < 2. This is explained on 
one hand by sensitivity limitations of the infrared instru- 
ments previously used, such as ISOCAM and ISOPHOT 
on-board ISO, and on the other hand by the SCUBA sur- 
veys being mostly sensitive to L > 1O 12 L galaxies which 
were much more numerous at z ~ 2-3. Current knowl- 
edge indicates that a large fraction of today's stars were 
probably born at 1 < z < 2 (e.g., Dickinson et al. 2003; 
Calura & Matteucci 2003), before the global activity of 
star formation began to decrease rapidly down to z ~ 
(Lilly et al. 1996). 

The recently commissioned Multiband Imaging Pho- 
tometer for Spitzer (hereafter MIPS) on-board Spitzer 
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provides a unique opportunity to address this issue. For 
sources at z > 1.5, its 24 (im band is particularly suit- 
able for detecting the possible redshifted emission from 
the 8 /xm broad-band feature commonly seen in infrared 
galaxy spectra, while the main 15 /xm filter on ISOCAM 
was rather more adapted for detecting this feature at 
z ~ 0.8. Here we assess the capabilities of the MIPS in- 
strument to shed light on a yet unexplored side of galaxy 
evolution, and present our preliminary results on the de- 
tection of 1 < z < 2.5 infrared galaxies at 24 /xm. This 
work is one facet of a larger set of publications also ad- 
dressing source number counts at 24, 70 and 160 /xm (Pa- 
povich et al. 2004; Dole et al. 2004), IR model interpre- 
tation (Lagache et al. 2004), as well as the relation be- 
tween the MlPS-selected objects and other populations 
such as the SCUBA/MAMBO and VLA galaxies (Egami 
et al. 2004; Ivison et al. 2004; Serjeant et al. 2004) and 
X-ray sources (Alonso-Herrero et al. 2004; Rigby et al. 
2004). Throughout the paper, a ACDM cosmology with 
H = 70 km s _1 Mpc _1 , Q m = 0.3 and ft A = 0.7 is as- 
sumed. 

2. OBSERVATIONS AND DATA ANALYSIS 

We performed 24 /xm observations 7 with MIPS (Ricke 
et al. 2004) in a 5'x5' area of the Lockman Hole (LH, a = 
10 h 51 m 56 s , <5 = 57°25'32", J2000) as well as in a 2°xl0' 
region of the Extended Groth Strip (EGS, a ~ U h 19 m , 
S ~ 52°48', J2000). Data in the Lockman Hole were ob- 
tained in the MIPS "photometry mode", with a series of 
dithered exposures allowing a high frame redundancy and 
efficient cosmic ray rejection (see Egami et al. 2004 for 
more details). Observations in the EGS were performed 
with the MIPS "slow scan" technique, an observing mode 
that allows coverage of large sky areas with high efficiency. 
Data reduction was carried out using the MIPS Data Anal- 
ysis Tool (Gordon et al. 2004). Effective integration times 
of ^250s and ^450s per sky pixel were reached respec- 
tively for the Lockman Hole and EGS data, leading to a 
24 /xm la rms -35/xJy (LH) and -25/xJy (EGS). The as- 
trometry of the Groth Strip final mosaic was refined using 
bright objects also listed in the 2MASS catalog (Jarrett 
et al. 2000), while source positions in our smaller coverage 
of the Lockman Hole were matched against an abundant 
set of ancillary data (see below) . We estimate the absolute 
astrometry accuracy to be better than 2". Point source ex- 
traction and photometry were performed as described by 
Papovichet al. (2004). ' 

Both fields were also observed with the Spitzer InfraRed 
Array Camera (hereafter IRAC, Fazio et al. 2004) at 3.6, 
4.5, 5.8 and 8/xm, down to ha rms ~1.2, 1.1, 6.5 and 
8.0 /xJy respectively. A detailed description of the IRAC 
data analysis is presented by Huang et al. (2004). In addi- 
tion, we retrieved H ST/WFPC2 images of the EGS from 
the public database of the Deep Extragalactic Evolution- 
ary Probe (DEEP) team 8 , and a set of ancillary observa- 
tions of the Lockman Hole was collected in the optical and 
the near-infrared (NIR) to supplement our Spitzer pro- 
gram. For the latter, we made use of archival i3-band 
images from the Wide Field Camera (WFC) of the Isaac 
Newton Telescope (INT), while VRIJHK S data were ob- 

7 Spitzer Prog. IDs 8 & 1077. 

8 see http://deep.ucolick.org 



taincd with the UH8K camera on the Canada-France- 
Hawaii Telescope (UT-bands), the Suprime-CAM instru- 
ment on the Subaru Telescope (i?-band) , TIFKAM on the 
2.1m Telescope of Kitt Peak (.ff-band) and Omega-Prime 
on the 3.5m Telescope of Calar Alto (JK s -ba,nd). Limit- 
ing Vega magnitudes (3cr level) of 26.0, 24.0, 27.0, 24.0, 
22.0, 20.5, 20.0 magarcsec -2 were respectively derived in 
the BVRIJHKg filters. Source detection and photometry 
were carried out with Sextractor (Bertin & Arnouts 1996). 
For the Lockman Hole data, the full-width at half max- 
imum of the PSF goes from ~1" in the optical/NIR bands 
to 2" and 6" in the IRAC and MIPS images, which gave 
us confidence in identifying the respective counterparts at 
each wavelength (see an illustration in Fig. 1). Merged cat- 
alogues among the different filters were then built, starting 
with matching the coordinates of the 24 /xm sources in the 
deep i?-band image. Following Kron (1980), we estimated 
the total flux for each object detected in the R image 
using an elliptical aperture derived from an isophotal fit 
around the given source, and translated those apertures to 
the other optical and NIR data to obtain integrated fluxes 
for the corresponding filter (in all cases the apertures were 
large enough to enclose the PSF profile). IRAC and MIPS 
fluxes were derived respectively with circular aperture pho- 
tometry and PSF fitting (see Huang et al. 2004; Papovich 
et al. 2004). 

3. RESULTS: INFRARED GALAXIES AT Z > 1 

3.1. Spectroscopic redshifts 

Our source position catalog in the EGS region was com- 
pared to the DEEP1 database of galaxies with spectro- 
scopic redshifts (Weiner et al., in preparation). Though 
most of objects from this survey clearly lie at < z < 1, 
we identified from our 24 /xm observations nine galaxies 
at 1 < z < 1.3, and two other sources at 0=1.43 and 
2=1.74 whose redshifts were respectively derived from Ha 
(Hopkins et al. 2000) and Fe II absorption lines (B. Weiner, 
priv. communication). From their observed 24 /xm fluxes 
and using the luminosity-dependent spectral energy dis- 
tribution (SED) templates from the model of Lagache 
et al. (2004), we then estimated their infrared lumi- 
nosities between 8 and 1000 /xm. Even though the ex- 
act relation between the mid- and total IR galaxy emis- 
sion is not yet well understood, there are some indica- 
tions that the mid- to far-infrared correlation still holds 
in the distant Universe (Appleton et al. 2004), and the 
L24/im/LiR ratio anyhow varies by less than a factor of 
~2-3 for most model SEDs at high reshifts (e.g., Pa- 
povich & Bell 2002). With this caution in mind, we found 
4 sources with 4xl0 11 L Q < Ltr < 10 12 L Q , 6 galaxies with 
10 12 L Q < Lir < 3xlO 12 L , and one hyper luminous sys- 
tem characterized by Lir=2.5x1O 13 L . Given that the 
objects from the DEEP1 database are optically bright, we 
note that these luminous sources likely have higher opti- 
cal/IR flux ratios than those of more obscured infrared 
galaxies. This underlies the importance of obtaining near- 
IR and mid-IR (e.g., IRS/Spitzer) spectroscopy as well as 
accurate photometric redshifts to study all types of IR- 
luminous objects. 
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3.2. Photometric redshifts 

For distant galaxies, our analysis of SEDs for EGS 
sources with spectroscopic redshifts indicates that the 
combination of the 4 IRAC bands with optical and NIR 
data allows locating with good accuracy the spectral fea- 
ture produced both by the stellar H~ opacity minimum 
at rest-frame 1.6 /xm and the global shape of the underly- 
ing stellar continuum. This characteristic of galaxy SEDs 
provides a key indication for photometric redshift deter- 
minations (Huang ct al. in preparation, see also Simpson 
& Eisenhardt 1999; Sawicki 2002). Figure 2a illustrates 
this result in the case of two representative sources of our 
sample, whose broad-band SEDs lead to photometric red- 
shifts in good agreement with those determined from op- 
tical spectroscopy. 

Relying on this new approach, we used our optical/NIR 
and IRAC data to locate 24 (im selected sources with 
photometric redshifts z > 1. The analysis of the com- 
plete MIPS data set in the EGS is still in progress. 
We present in this preliminary study the results we ob- 
tained for the Lockman Hole. Over the 5'x5' area cov- 
ered by MIPS, we identify at least 25 objects charac- 
terized by a rather prominent stellar bump revealing an 
unambiguous redshift larger than 1 . Figure 2b shows the 
BVRIJHK /IRAC/MIPS broad-band SEDs of several of 
these sources spanning the 1 < z < 2.5 redshift range, with 
the best galaxy template fitted to our data. These tem- 
plates were built by combining the optical/NIR spectral 
range emission modeled by Devriendt et al. (1999) with 
the infrared SEDs from the model of Lagache et al. (2004). 
For several objects showing a rather high MIPS to IRAC 
flux ratio, an additional contribution from a hot dust con- 
tinuum (i.e., Very Small Grains) was arbitrarily added (see 
Gallais et al. 2004 for more details on the method) . In this 
first attempt to locate z > 1 galaxies using the 1.6 /im stel- 
lar bump, we relied on a simple "by-eye" fit which clearly 
gives a sufficient constraint (i.e., Az~0.2) for the cases 
shown in Fig. 2. We emphasize that the few SEDs shown 
in Fig. 2 were not chosen to be globally representative of 
our galaxy sample at similar redshifts, but correspond to 
24 /im sources detected at least in the -R/IRAC bands and 
for which a photometric redshift could be reliably obtained 
from the spectral feature produced at rest-frame 1.6 /im. 
This naturally biases the selection against sources with 
a rather shallow stellar bump like the active nuclei (see 
Alonso-Herrero et al. 2004). This issue will obviously have 
to be taken into account in future studies of more complete 
samples. 

Table 1 summarizes general properties of our sample in 
the Lockman Hole. As in the EGS, the estimated red- 
shifts and observed 24 /im fluxes indicate infrared lumi- 
nosities characterized by Lir > 10 11 L Q . The rest-frame 
optical and NIR properties reveal rather red colors (i.e., 
observed R — K > 4) originating cither from an evolved 
stellar population or a young reddened starburst. They 
also indicate high luminosities in the rest-frame near- 
infrared, which suggests the presence of rather massive 
objects (M >M*). 

Finally, we note that only one of these sources is de- 
tected in the deep X-ray observations of the Lockman Hole 
by XMM (Hasinger et al. 2001). In conjunction with the 
galaxy-like templates fitting their optical SEDs, the X-ray 



non-detections suggest that these galaxies are mostly dom- 
inated by star-forming activity, even if the presence of an 
active nucleus can not be fully excluded by the data (see 
Alonso-Herrero et al. 2004). 

4. DISCUSSION AND CONCLUSION 

We stress that our reported number of 24 /xm-selected in- 
frared galaxies at z > 1 is obviously a lower limit on their 
density at these redshifts, since we only selected sources 
characterized by a well-defined stellar bump in the rest- 
frame near-infrared. This limit agrees with the predictions 
of the various models dealing with infrared galaxy evolu- 
tion (e.g., Lagache et al. 2004). 

Our selected sources exhibit a rather wide range of MIPS 
to IRAC flux ratios and optical/NIR SED shapes. This 
suggests a possibly large diversity in the properties of in- 
frared galaxies at high redshift. It is also interesting to 
note that the most distant sources of our sample tend to 
have extremely red colors from the B to the 24 /im band, 
with a global SED best fitted by ULIRG templates such as 
the one from Arp 220. Finally, our data suggest the pres- 
ence at high redshift of the broad-band emission features 
often seen in the mid-infrared spectra of local galaxies and 
commonly attributed to large molecules transiently heated 
by UV photons. Even though we can not completely ex- 
clude a steeply rising continuum accounting for the 24 /im 
emission, the observed MIPS/IRAC flux ratios are indeed 
difficult to reproduce without a significant contribution 
from these features (see Fig. 2). A similar conclusion is 
reached by Lagache et al. (2004) based on the modeling of 
the MIPS 24/70/160 /im source number counts. 

Since the late 1990's, there have been a wealth of discov- 
eries at very high redshift (e.g., z > 3), while little progress 
has been made in the 1 < z < 2 redshift range (but see 
Stcidel et al. 2004). Yet, galaxy evolution at 1 < z < 2 
played a central role in the cosmic growth of structures 
and the assembly of galaxies (Ellis 2001). MIPS provides 
a unique opportunity to explore the star formation and 
black hole accretion history at these redshifts. Compared 
to previous models for instance, the new predictions by 
Lagache et al. (2004) based on the MIPS source number 
counts indicate a high density of z<~1.5 LIRGs/ULIRGs, 
the first of which have just been identified in this paper. 
The MIPS deep surveys will reveal the evolution of lu- 
minous and massive galaxies at 1 > z > 3, bridging the 
gap between the different populations revealed by ISO 
and SCUBA. Combining with sources at similar redshifts 
but selected with different techniques (e.g., Balmer-break 
galaxies) will also give us new clues on galaxy evolution. 
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Fig. 1. — An illustration of our band-merging process, showing a source detected with MIPS at 24 (im (a) and its counterparts at 8/im (b), 
3.6 jim (c) and 0.6 ^m in the R— band (d). Contours in these 3 last panels correspond to the 2Afj,va emission. The centering accuracy at this 
wavelength (< 2"), combined with the sharp PSF of the IRAC data between 3 fim and 9 (im allow unambiguous identifications. The four 
images have a similar field of view, with the scale indicated in the top-left panel. We estimate for this source a photometric redshift z~1.8 
(see Sect. 3). 
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Fig. 2. — Optical/NIR/IRAC/MIPS broad-band photometry ('*' symbols) of 24fim selected sources with spectroscopic (a) and photometric 
(b) redshits z > 1. In each case, the dashed curve represents the best redshifted SED template fitting the data (see Sect. 3.2 for more details). 
The corresponding photometric redshift and the source coordinates (J2000 epoch) arc mentioned to the right of the panel. Flux uncertainties 
are represented with vertical solid lines. A 2=0 starburst-likc spectrum is added at the bottom for reference. Note how the rest-frame 1.6 fJ.ni 
stellar bump gives a clear constraint on the redshift determination (JJ. symbol). 
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Table 1 
General properties of our selected sources in the Lockman Hole 



Rcdshift 


type™ 


N 6 


Lir( 


:io n L r 


M K c ' d 


1.0- 1.5 


Red Sp. 


>3 




4! 


nt o-24.8 
zo - -25.6 




Starb. 


>H 




6j° 


-24.8lgg;? 




VLIRG 


>3 




18?? 


9 c 7 -24.9 
z0 -'-26.3 


1.5 - 2.0 


Red Sp. 


>1 




1 


-26.4 




VLIRG 


>7 




26 f 


_ 254 -24.6 
zo.<i_ 27 


2.0 - 2.5 


VLIRG 


>3 




5818 


-26.0I^;4 



a Source classification — Red Sp.: IR-quicscent early- type spi- 
ral galaxy with a red optical/NIR color. Starb.: M82-likc IR 
starburst associated with a late-type spiral. VLIRG: Very Lu- 
minous IR Galaxy characterized by a high rcst-framc IR/NIR 
ratio and a rather elusive 4000A break. 

b Minimum number of selected sources per category. 

c assuming Ho = 70 km s Mpc -1 , f2 m =0.3 and £l\ = 0.7. 
The mean value is given, along with the range of measurements. 

d Rest-frame A"-band absolute magnitude. 



